The Widom-Dyson constant for the gap probability in 

random matrix theory 



P. Deift 

Courant Institute of Mathematical Sciences, 
New York, NY 10003, USA 

A. Its 

Department of Mathematical Sciences, 
Indiana University - Purdue University Indianapolis 
Indianapolis, IN 46202-3216, USA 

I. Krasovsky 
Department of Mathematical Sciences, 
Indiana University - Purdue University Indianapolis 
Indianapolis, IN 46202-3216, USA 
and 

Department of Mathematical Sciences 
Brunei University 
Uxbridge UB83PH 
United Kingdom 

X. Zhou 

Department of Mathematics, 
Duke University, 
Durham, NC 27708-0320 



1 



1 Introduction 



In this paper we consider an asymptotic question in the theory of the Gaussian Unitary 
Ensemble of random matrices ~JJ. In the bulk scaling limit, the probability that there are no 
eigenvalues in the interval (0,2s) is given by P s = det(J — K a ), where K s is the trace-class 
operator with kernel 

sin(x - y) 
ir(x - y) 

acting on L 2 (0, 2s). We are interested particularly in the behavior of P s as s — > oo. 
In 1973, des Cloizeaux and Mehta j2j showed that as s —>■ oo 

lnP s = -£-~lns + c + (l), (1) 

for some constant c. In 1976, Dyson '3J showed that P s in fact has a full asymptotic expansion 
of the form 

s 2 1 a 1 a 2 

7-4 lns + Co + 7 + ^ 

Dyson identified all the constants Cq, ai, a 2 , ■ ■ ■■ Of particular interest is the constant c , 
which he found using earlier work of Widom (see jl] and below) to be 



lnP s = - — -- In s + c + — + -£ + ■••■ ( 2 ) 



c = ^ln2 + 3C / (-l), (3) 

where ((z) is the Riemann zeta-function. 

The results in [2] and [3] were not fully rigorous. In j^j, Widom gave the first rigorous 
proof of the leading asymptotics in in the form 

InP. = -£(1 + 0(1)). (4) 

In subsequent work [HI E] , which also included the multi- interval generalization, the form (j2J) 
of the full asymptotic expansion was verified rigorously, together with the correct constants 
ai, a 2 , .... The expression (JSJ) for the constant Cq, however, remained unproven. This was 
because the methods in jSJ E] and [7] naturally computed the asymptotics of (d/ds) In P s , 
and the constant of integration remained undetermined. 

Recently, two proofs of Q were given independently in the literature in jH] and [9*1 ITU]. 
The methods in the papers jSj and [TO] are very different. Our goal in this paper is 
to give a third proof of (j!JJ), which is closely related to the proof in jS], but as explained 
below, does not rely on certain a priori information. This means that our approach has 
the potential advantage of being applicable to other problems involving the computation of 
critical constants, where a priori information may not be available (see, e.g., |T""]). 
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One way that one might try to evaluate Co is to express 



lnP s = lndet(/-fQ = [ ^-tr ln(J - r)K s )dr) = - [ tr ((I - nK^KMri (5) 

Jo d V Jo 

and then evaluate tr ((/ — r]K s )^ 1 K s ) asymptotically as s — > oo for each fixed 77 G (0, 1) using 
steepest descent methods as in jTHj, for example. However, it turns out that the asymptotics 
of tr ((/ — 7]K S )~ 1 K S ) as s — > 00 have a different form for 77 < 1 and rj = 1. This means that 
one must integrate the asymptotics in (jSJ) over a boundary layer as 77 —>■ 1, a difficult task 
which we have so far been unable to perform. On the other hand, for < 7 < 1, we can 
indeed use (J5J) in the form 



lndet(J - jK.) = - fti ((/ - i]KsT l K t 
Jo 



together with the Riemann-Hilbert /steepest-descent method to compute the asymptotics of 
lndet(J — jK s ) as s — > 00, so reproducing the results in fH\ IT2]. 

As mentioned above, Dyson's computation of Co in j^j is based on an earlier calculation 
of Widom In [I], Widom considered, in particular, the Toeplitz determinant D n (a) with 
symbol given by the characteristic function of the interval (a,2ir — a), < a < 71. Thus 
D n (a) = det(M i _ i )"7i > where M k = J^~ a e- lk9 d9/27T, k G Z. Widom showed that for a 
fixed a as n — » 00, 

\nD n (a) = n 2 lncos — - - In (^nsin — J + c + o(l), (6) 
where Cq is the constant (jHJ). What Dyson noted was that for a fixed s > 0, 

lim D n (-) = det(J - K s ) = P s (7) 



n^oo \ 77, 



and hence, if the error term o(l) in (jUJ) was uniform asn-> 00, « — > 0, an — > 00, one could 
conclude from (|6I7|) that Co in (0) is indeed given by 0- The main technical result in this 
paper, as in jS], is the proof that the error term o(l) is of the form 0(l/(nsin |)), which 
gives the desired uniformity. 

Whereas P s is the gap probability for the Gaussian Unitary Ensemble in the bulk scaling 
limit, we note that D n (a) is the gap probability for the Circular Unitary Ensemble pQ. 
Formula (|IJ) is the scaling limit for this probability, and the fact that the limit also gives P s 
is a well-known universality property. 

In jH] , the author uses steepest descent methods to show that for e > fixed, there exists 
a (large) positive constant so such that 

d , _ . . n 2 a 1 a ^ ( 1 \ 

— In DJ a) = tan cot - + 9/ , \ (8) 

da K J 2 2 8 2 V^sin 2 (a/2)y v ; 
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for all n > so and ^ < a < it — e. Integrating (jHJ) over (a, ao), 2so/n < a < ao < it — e, 
one obtains 



,2, „„ « 1 ,„^ a „2i «0 



ln_D n (a) = ln_D ra (ao) + n In cos — — - In sin — — n In cos— + 

1 1 ■ a o , /n 
- In sm h C 



(9) 



4 2 V^ sin ("/ 2 )/ ' 

Using Widom's result © for fixed ao, one obtains for — < a < tt — e, n > sq, 



\nD n (a) = n 2 lncos^ - ~ In (rising) + c + O ( — .) ) + 8 n , 
2 4V 2/ \nsin(Q;/2)/ 



(10) 



where S n — > as n — > oo. For any fixed s > So, one sets a = 2s/n, and then using ((7|) and 
letting n — > oo, one obtains 

lnP s = -y-iln S + c + oQ, (11) 

which proves ©• 

In this paper, we will derive an improved version of Q, viz., 

\nD n (a) = n 2 lncos^ - \ In fnsin + c + O ( — — ) , (12) 
2 4V 2/ \n sm(a/2)/ 

for ^ < a < n — e, n > sq, where Sq is again a (large) positive constant. Our proof of (|12j) 
is direct and does not rely on Widom's result ©• The proof is based on the following two 
principles: 

(i) Asymptotics of D n (a) as a —>■ ix and n is fixed. 

(ii) Asymptotics of the solution of a regularized version of the [7] Riemann-Hilbert problem 
(see below) uniform for 2s /n < a < ir. 

The solution of problem (i) is based in turn on the analysis of the standard multiple- 
integral representation for D n (a). The solution of problem (ii) is based on a mapping of 
the original Riemann-Hilbert problem posed on the arc —a < 9 < a of the unit circle 
to a problem on the fixed interval [—1,1]. The analysis then proceeds via the steepest 
descent method for Riemann-Hilbert problems introduced by Deift and Zhou in [THj and 
further developed in [THj, (2H|, and also in [21]. This gives an asymptotic expression for the 
logarithmic derivative (d 2 /da 2 ) In D n (a). Formula ()21|) below together with its integrated 
version (jlHHj) . plays a key role in this paper. 

Note that in contrast to [5] and [7], where the analysis of the derivative (d/ds) \nP s fails 
to identify the constant Co, we may now integrate (d 2 /da 2 ) In D n (a) from a — > it, and the 
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limit at tx is determined by step (i). The result is the expression ([T2]). By contrast, in [3] 
and [7] there is no convenient point So from which we can integrate and then use to extract 
the relevant asymptotics. The key device that makes our method work is the $-RH: In 
particular, we note that the 11-element in the jump matrix for the $-RH (see fJHljl et seq.) 
is uniformly small as n — > oo, for all 2so/n < a < tt, sq » 1, and for all A in a compact 
subset of (—1, 1). It is this uniformity in a as n — > oo that makes it possible to control the 
integration from a = n to a = 2so/n. 

In Section 2 we analyze step (i), and in section 3, step (ii). Finally, in Section 4, we 
prove JX2J). 

2 Step (i). Multiple integral analysis. 

For the analysis of D n (a) as a — > tc we use the multiple integral (e.g., [TTf fTHj) 

where C a is the arc a < 9 < 2tt — a of the unit circle. The integrals are taken from a to 
2n — a. Setting 

a — 7r - (3, (3 > 0, 

and 

9j = ii + [3xj, 

we rewrite f!13|) as follows: 

D n { a ) = ^L^fr £••• f J] \e i ^-e i ^dx 1 ...d Xn . (14) 

Observe that 



1 l<j<k<n 



\ [ \e ifke ' -e^f = (3< n ^[ \ Xj - x k \ 2 + O n {(3 2 ) 

l<j<k<n \l<j<k<n 



Hence we arrive at the relation 



jDn(a) = T^W^" 2 1 / " / II I'r ^l 2 ^ 1 • • • dXn + °"^ 2 ) f • ( 15 ) 

I ) n - [7-1 l<j<k<n J 

The multiple integral in this formula can be expressed in terms of the norms of the Legendre 
polynomials. Indeed (see, e.g., [T7]) 

/i „i n-1 

" / I f l x i _ ^fc| 2 ^i • • • dx n = ii\ Y\ h k , (16) 

1 l<j<k<n k=0 
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where h n are the normalization constants of the monic polynomials orthogonal on the interval 
[—1,1] with the unit weight: 

p n (x) = x n + . . . , J p n (x)p m (x)dx = h n 8 nm . 

Let P n {x) denote the standard Legendre polynomials [T7j. Since 

(2*0! n , r^u„_ 2 



Pn(x) 



-x n + ... and j P%(x)dx 



2 n (n\) 2 J_, ny ' 2n + l 

we conclude that 

2 n M) 2 

Pn{x) = • P n (x) 



and 



1 2 , . , 2 2 >!) 4 /- 1 _ a . . . 2 2 >!) 4 2 



[(2n)!] 2 y_ x " v y [(2n)!] 2 2n + l' 
This leads us to the following representation of D n (a) in the neighborhood of a = 7r: 

ln£> n (a) = n 2 ln/?-nln2vr + lnA n + O n (/? 2 ), a = 17-/3, (17) 



where 

a = Y[ h , = Y[ 

\(2kW 2k + 1 

For later reference, note that the asymptotic relation (|T7j) is clearly different iable, for fixed 
n, with respect to a. Also, for fixed n, the term O n (f3 2 ) — > as /3 = 7r — a — > 0; no claim is 
made here about the behavior of O n (/3 2 ) as n — > oo. 

Widom's constant, Co = ^ In 2 + 3('(—l), is generated by the quantity A n . In fact, it is 
shown in j4j, using results from classical analysis, that 

A n = e co n- 1/4 (2rr) n 2- n \l + o(l)), n -> oo. (19) 

The appearance of the zeta function is due to the presence of the products of factorials. 
Indeed, 

n n n 

In fc! = V] In fc! = (n + 1) In n! - V] fc In fc, 
fc=i fc=i fe=i 
and one can expect that the asymptotics of the sum on the r.h.s. of the last equation is 
related to £'(— 1). The exact relation (see again 0j) reads as follows: 

n fl 1 1 \ 11 

^ In fc = I -n 2 + -n + — J In n - -n 2 + — - C'(-l) + o(l). 

fc=i ^ ' 

Applying this formula and the asymptotics of the Gamma- function to (|T8j) yields (fT9"j) . 
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Figure 1: Contour for the m-RH problem. 



3 Step 2. Riemann-Hilbert analysis. 

Denote the complement of C a in the unit circle by Y a = {—a < 6 < a} traversed coun- 
terclockwise (see Figure 1). Let m(z) = m(z;n,a) be the solution of the following 2x2 
Riemann-Hilbert problem posed on r a : 

• m(z) is holomorphic for all z ^ T a 

• m(oo) = I 



Here, as usual, m+(z) (respectively, m-(z)) are the L 2 boundary values of m(z') as z 1 — > 
z E T a non-tangentially from the "+" side {\z\ < 1} (respectively, "— " side {\z\ > 1}). We 
shall refer to this Riemann-Hilbert problem as the "m-RH problem" . 

Theorem 1 Let < a < it and n > 0. Then the m-RH problem has a (unique) solu- 
tion m(z; n, a), and the Toeplitz determinant D n (a) is related to m(z; n, a) by the following 
differential and difference identities: 




a 



D n (a) 



mn(0; n, a), 



(20) 



da 2 



In D n (a) 



[m 12 (0;n, a)] 2 . 



(21) 
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Remark As a function of z, m(z; n, a) has a continuous extension up to the boundary 
r a , apart from the two end points e ia and e~ ia , where it has logarithmic singularities. 
Moreover, m±(z) admit analytic continuations into a neighborhood of every point z of the 
open arc r Q = T a \ {e m , e~ ta }. Note also that detm(z] n, a) = 1 by a standard calculation. 
These properties of m(z; n, a) are inherited by solutions of the transformed Riemann-Hilbert 
problems introduced below. 

Theorem 1 was proved in [7] 1 (cf. Eqs (6.14) and (6.82)) using standard techniques from 
the theory of integrable systems: derivation of the relevant Lax pair, identification of D n (a) 
as the relevant tau-function etc. The differential identity (|21|) will be of central importance 
for the analysis below. 

A standard calculation shows that the m-RH problem has no solution for a = it. How- 
ever, as we now demonstrate, the m-RH problem can be regularized for all a in the range, 
including a = ir, by a simple sequence of transformations. 



3.1 Mapping onto a fixed interval. 

For < a < tt, the linear-fractional transformation, 

a z—1 1 + iAtanf 

A = -zcot , z = (22) 

2 z + V 1-iAtanf v ; 

maps the arc T a onto the interval (—1, 1) and transforms the m-RH problem to the following 
Riemann-Hilbert problem posed on the interval (—1, 1) traversed from —1 to 1 (see Figure 
2): 

• Y(X) is holomorphic for all A ^ [—1, 1] 

• Y(oo) = I 

1+tA tan 



>'-(A. = V + (A)| , ^T**' 1 Ae(-I.l, 

1-iAtanfy U 

We shall refer to this Riemann-Hilbert problem as the "F-RH problem". The relation 
between the F-RH problem and the original m-RH problem is given by the equation 

m(z; n, a) = F _1 (— i cot — ; n, a)Y(X(z); n, a). (23) 



1 There are some differences from the notation in 0, namely, our contour is r a instead of C {T a rotated 
by 7r) in and we make the following choice for the functions fi, gi which build up the kernel: j\ = z n / 2 , 
H = z- n '\ ,9i = z-"/ 2 /(2™), g 2 = -W 2 /(2™). 



® 



^ + 

-1 " 1 

Figure 2: Contour for the F-RH problem. 

The F-RH problem is still irregular at a = it. Indeed, the function 0(A, a) = fzixlwa/^ 
is discontinuous at (A, a) = (0, tc). We have for the jump matrix Jy(A, a) of Y(X) as a — > 7r: 



Jy(A,7T) 








1,1 A ^ 



A = 0, 



which demonstrates the difficulty for odd n. For even n, however, the jump matrix Jy(A, tc) 
is continuous and constant throughout the whole interval (—1,1). This implies the solvability 
of the y~RH problem at a = it for even n; in fact, one easily checks that 

y(A;»=2»,»)=(; i) (i ^>^)(; iy. 

However, regardless of the parity of n, the convergence of Jy(A, a) to Jy(A, 7r) is not uniform 
in A, and this creates a significant difficulty in the direct analysis of the behavior of the 
solution Y (A; n, a) near a = it. 

As we now show (see the $-RH problem below), the ^-problem can be regularized by 
performing one more step which is familiar in the formalism of the nonlinear steepest descent 
method. 



3.2 (/-function transformation. 

Following the nonlinear steepest descent method for Riemann-Hilbert problems (see, e.g. 
[21]), we introduce the following "g" function: 



tW s i+ ;^!f f . (24) 



1 + iX tan | 



This is essentially the ^-function of section 6 of [7j written in the variable A (see equation 
below). It possesses the following characteristic properties: 



(a) g(X) is holomorphic for all A ^ [—1,1]. Here we fix the square root by the condition 

VX 2 — 1 ~ A, A — ► oo. 
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(b) g{X) ^ for all A ^ [—1,1]. At the points A = — i cot | (or z = oo) and A = i cot f (or 

z = 0) the values of the function g(X) are: 

,. a a 2 a 

q{— i cot— ) = 1 and g(z cot — ) = cos — = k. 

(c) The boundary values of g±(X), A G [—1, 1] satisfy the following equations: 



9+9- = « 



1 — iA tan |- 
1 + iA tan § 



and 



9+ 



sin 



sin 



(d) The behavior of g(X) as A — > cxo is described by the asymptotic relation 

^( A ) = cos| + oQ 



(25) 



(26) 



(27) 



(2? 



It is worth noticing that 



z + 1 + - e ia )(^ - e~ ia ) 



2z 



(p(z). 



(29) 



If one changes 1 to —1 in the numerator, then (p(z) becomes the ^-function of section 6 of 
[7j. The change of sign is due to the fact that the Riemann-Hilbert problem considered in 
[7j is defined on the arc C = e m T a rather than on T a (cf. footnote 1 above). 

Equation (J27j) has an important consequence. Fix < 5 < 1 and < a < n. Then the 
following inequality holds: 



9+ 



9- 



<e <l, AG [-1 + 5,1-5], 



(30) 



for some e$ = £o(5, a) > 0. Of course, for all A G (—1, 1) and a G (0, 7r), we have 

9+ 



9- 



< 1. 



Following the steepest descent method, we transform the original Riemann-Hilbert prob- 
lem by the formula 

(31) 



where a 3 







Y(X) ^ $(A) = Y(X)g- na3 K^ 3 



is the third Pauli matrix. From the properties of the ^-function listed 



-1 

above, it follows that the matrix function $(A) = $(A; n, a) is the solution of the following 
Riemann-Hilbert problem, which we shall refer to as the "$-RH problem" : 
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$(A) is holomorphic for all A ^ [— 1, 1] 



$_(A) = $+(A) 



-1 



A G ( — 1, 1; 



In view of ()23|). the original function m(z) is related to the solution $(A) by the formulae: 



mlz: ^, a 



K2 



cot -; n, a)$(A(z);n, a)c/ nff3 (A(z)K 



K2 ff3 $ *(— jcot ^;n, a)<&(A(z);ra, a)v9 n<J3 (z)/t ^ 



(32) 



As indicated earlier, the $-RH problem is regularized. Indeed, note first that the jump 
matrix for the $-RH problem is now continuous for all A G [—1, 1] and a G [0, 7r] with the 
end point a = tt included. Moreover, for all < a < tt, the $-RH problem is (uniquely) 
L 2 -solvable, by the following argument. Consider the $-RH problem as defined on the whole 
real line with discontinuities at ±1 (the jump matrix outside (—1, 1) is the identity). The 
limiting value of the jump matrix at these two points from inside the interval (—1, 1) is 
2 — 1 \ 

whose only eigenvalue is 1 ^ (— oo, 0]. By Theorem 5.16 of [2H], the RH problem 



1 

is L 2 Fredholm. Now by (the proof of) Theorem 9.3 in [26J, a Fredholm RH problem with 
a jump matrix v on K is L 2 -solvable if v + v* > everywhere, and v + v * > on a set of 
positive Lebesgue measure. These conditions are clearly satisfied in our case. Therefore the 
$-RH problem is L 2 -solvable. 

Theorem [T] and equation (|3*2*|) yield representations of the Toeplitz determinant D n (a) 
in terms of the solution of the $-RH problem: 



D 



n+l 



a) 



dot' 



D n (a) 
In DJa) 



Q(n, a) cos 



2d 



a 



n 



sin a 



■2 _ A ( n >«)> 



where 



and 



0(n, a) 
A(n, a) 



(33) 
(34) 



01 Oi 

i cot — ; n, a)$(z cot — ; n, a) 

2 2 j ii 



01 OC "I ^ 

i cot — ; n, ai)Q>(i cot — ; n, a) 

2 2 J 12 



(35) 
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3.3 Asymptotic analysis of the <3>-RH problem. 



By standard arguments, using inequality (|30|) . one expects that $(A) is approximated by the 
function 

" / 3(A)+/3- 1 (A) /3(A)-/3~ 1 (A) 

N(X)= I 2 * ] , P(X) = j] , (36) 



PM-p-ijX) /3(A)+/3- 1 (A) / ' ' \A + 1 



2i 2 

/S(oo) = 1, 

which solves the model Riemann-Hilbert problem: 

• iV(A) is holomorphic for all A ^ [—1, 1] 

• JV(oo) = / 

. JV.(A) = JV + (A) I I, A €(-1,1) 

In order to estimate the precision of this approximation, we need to consider the $-RH 
problem for A near ±1. The following result, which allows for complex values of a in a 
neighborhood of a = tt, is basic for our analysis (the need for this complex extension will be 
apparent towards the end of the paper, see (jl26J) below). 

Theorem 2 Let 5 be a positive number less than 1/4. Introduce the domain 



ftW =C\(UUU), 

were U (U) denotes the open disk of radius 5 centered at 1 (respectively, —1). Let also e 
be a positive number less than n — 2 and denote T> £ (tt) the disk in the a-plane of radius e 
centered at a = n. Set 



p = n 



a 

sin — 

2 



Then, for 5 and e sufficiently small, there exists s > such that for all a G [— , 7r — e\ U 
V e (7i), and n > so, the solution of the ^-RH problem exists (and is unique) and satisfies the 
estimate 

$(A;n,a)= (l + O ( * Y) N(X), p - oo, (37) 

uniformly for A G f^ 2 ^ and a G [^,7f — s] UZ> e (7r). Moreover, this estimate can be extended 
to a full asymptotic series in inverse powers of p; in particular, the order p~ 3 extension of 
PTp reads: 

$(A; n, «) = (/ + Ri(X) + R 2 (X) + R r (X)) N(X), (38) 
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where 



16ip 



R 2 (X) 



2 8 p 2 



.- A \ / 

_ ( 1 8i 
A V -8i 1 



i? r (A) = O 



1 A l 
i + TTa 



AgC\((/U[/), (39) 



+ 



1 + A 



OO, 



AgC\([/U[/), (40) 



(41) 



:i + iadp 3 

uniformly for A G fi (25) and a G [^a, 7r - e] U V £ (n). 

Remark 1 The last statement (|41|) means that there exist positive constants C and sq 
depending on e and 5 only, such that 

|i2r(A)| < 



c 



(42) 



VA G fi (25) , Va G 



2s 



n 



-, 7T — £ 



U 2? s (7r), Vn : s < ri. 



We shall also assume that e is small enough for the inequality, 



a 



sm 



> 



to take place for all a G T> £ (ir), and hence 



1 



P > 2 S o, 



(43) 



(44) 



n. 



for all a G [^f, vr - e] U V £ (ir) and s < 

Remark 2 Part of the assertion of Theorem 2 is that the solution of the $-RH problem 
exists and is unique for all a G [2so/n, n — e] U V E (ir) and n > so with sq sufficiently large. 
This is all we need in the analysis that follows; however, the solution of the $-RH problem 
actually exists and is unique for all a G [0, it — e] U T> e (Tc) and all n > for some (possibly 
smaller) e > 0. Indeed, by the discussion following (|32p. the $-RH problem is solvable for 
all a G [0, 7r], n > 0, and also for all a G ZV(7r), < n < s for some e' > by continuity of 
the jump matrix at a = n. By Theorem 2, the $-RH problem is solvable for all a G T> e (jr), 
n > sq. Thus the $-RH problem is solvable for all n > on [0,7r — £i] U P £1 (7r), where 



mm e, e 



Remark 3 The local analyticity of the jump matrix of the $-RH problem implies that 
both boundary values of the function $(A) on (—1, 1), i.e. the functions <&±(x), admit the 
analytic continuation in the neighborhood of every point of the interval (—1 + 5,1 — 5). 

Proof of Theorem [21 We shall now construct parametrices in U and U which are 
solutions of the $-RH problem in these neighborhoods with the condition at infinity replaced 
by the requirement that they match N(X) at the disks' boundaries to leading order (cf. jZllHl)- 
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+ 





1 1 

Figure 3: Contours for the Y- and $-RH problems and parametrices. 



Consider the function 



/(A) 



1 + ?(A 2 - l) 1 / 2 sin(q/2) 
l-i(X 2 - l) 1 /2 s in(Q;/2) 



(45) 



which is analytic and has no zeros in U \ (1 — 5, 1] and U \ [—1, —1 + 5). (Note, however, 
that it is singular outside of these disks at A = — i cot(a/2).) The branch of the root is taken 
such that (A 2 — l) 1 / 2 > for A > 1. The function ()45|) has the following boundary values on 
(-1.1): 



/ +M = ;-^zg 3 ^ , /_(«)= / +( *)-', xe(-m). 

1 + Vl-r sm(a/2) 



(46) 



Consider first the neighborhood U. We look for a parametrix, an analytic function in U \ 
(1 — 5, 1], satisfying the jump condition of the $-RH problem on (1 — 5, 1), of the form 



P(A) = £(A)P(A)/(A)-^ 2 , A 6 U \ (1 - 8, 1], 



(47) 



where E(\) is a non-zero analytic matrix-valued function in U (which therefore does not 
affect the jump condition) to be chosen below so that P matches iV to leading order on the 
boundary dll. 

It is easy to verify using ()46|) that for P to satisfy the jump condition for the $-RH 
problem across (1 — 5, 1), P must satisfy the jump relation 



p.(x) = P-(x) 



1 

-1 2 



x e (1-5,1). 



(48) 



An appropriate matrix function satisfying this jump relation was constructed in [2j (cf. [Zj 
(4.79), (4.871)). 

For A G U \ (1 — 5, 1], define the analytic function 



«(A) = -ln/(A) 



(49) 
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= i y ( _i)*_L_ sin 2fc+i ^ (A 2 _ i)w/j. 

^ v ; 2fc + 1 2 V ; 

k=0 



Note that 



(A) = zv / 2sin-(A-l) 1/2 G(A), 
2 



UJ 



where G(X) is analytic in all of U, and 



1 2 



a 



G (A) = l + (A-l)( i --sin 2 -)+0((A-l 



for A near 1. Thus, 



a w(A) 



/(A) , Aetf\(l-*,1], 
u;(a;) + = e i7r u;(x)_, a; G (1-5,1). 

Furthermore, the function cu 2 (A) is analytic in all of U and 



a 



u 4 



a 



= 2ue m sin 2 - ( 1 + - - -u sin 2 - + 0(u 



u = A- 1. 



(50) 



(51) 



(52) 
(53) 



(54) 



The term 0(w 2 ) in ()51J54j) is uniform for all < a < it. In fact, the estimate (|54|) is uniform 
for a belonging to any compact set in the complex a-plane. Let us choose < e < n. Then, 
for sufficiently small 5, the asymptotic relation implies that 



(A)| > V5 



. a 
sin — 
2 



VXEdU, Va G [0,tt -e] UX> £ (tt). 



(55) 



Here T> £ (tt) is the disk in the a-plane of radius e centered at a = tt. 
Introduce the new variable 

( = e""nV(A). 
Note that the mapping A —>■ £ of U is one-to-one. 

^From ()54|) and ()55|) it follows that for S and e sufficiently small, the following inequalities 



(56) 



hold: 



and 



37T r- 37T 

T <ar gv /C< T , 



a 

sin — 

2 



(57) 
(58) 



VA G dU, Wa G [0, tt - e] U £> e (7r) 
Inequality (J5%|) together with imply the estimate 

iVci>^o>i, 



(59) 
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VA e dU, Va e 



2s 



n 



IT — £ 



UV £ (n), — <So <n. 



A function P(X) analytic in U\(l—5, 1] and satisfying ()48|) is given by the following expression 
in terms of Hankel functions (cf. [7j) where i/C = e~ l7T ^ 2 ncu(X): 



P(A) 



(60) 



Inequality (J57)l and estimate (JoT?|) allow us to use the standard expansion for Bessel 
functions and obtain the following asymptotics on the boundary dU: 



P(X) 



2 A-os/4 f 1 1 
7T \ / -/ 



' + ^(-2 -l)+#c(-4 (61) 



e nc;(A)(T3 e -i(7r/4)fT3 



A G <9£7, 



where the remainder P r (X) satisfies the uniform estimate 

C 



|P r (A)| < 



VA G dU, Va G 



2s 



n 



-, 71 — s 



| C |3/2> 



UV £ (n), — <So < n . 



(62) 



Here Cq is a numerical positive constant which comes from the universal asymptotic expan- 
sion of the Hankel function H^\y/() for \(\ > 1 and — 3tt/4 < arg y/( < 3tt/4. 

Now let us choose E(X) so that P matches N on dU to leading order in p, i.e., PN^ 1 ~ /. 
Clearly, we should take 



E(X) = N(X)e 



i(ir/4)cr 3 _ 



1/1 -i 



2 V 1 



<7 3 /4 



(63) 



It is easy to verify that E(X) has no branch point or singularity at A = 1. Hence E(X) is 
analytic in U. 

Thus, the parametrix in the neighborhood U is given by the expression: 



P(A) = N(X)e i{n/4)a3 



7T / 1 



2 3 / 2 V 1 * 



-I \ XT3/4, 



C CT3/4 x 



(64) 



where ( and /(A) are defined by (|56|49|45J) . 



/(A) 



-0-371/2 
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The construction of a parametrix in the neighborhood U is similar. In this case, instead 
of (HJ) we set 

S/(A) = ~ln/(A), (65) 
which is analytic in U \ [— 1, — 1 + 5). Thus 

e s{x) = /(A) i/2 ) AG ^\[_ lj _ 1 + ff ) j (66) 
= e~ i7r £(x)_, are (-1,-1 + 5). (67) 

We find the same power series expansion for ui 2 (X) as (|54|) with A — 1 replaced by —A — 1. 
We define the ( variable for A in U again by the equation 

( = e- in n 2 u 2 (X). (68) 

Note that for both the images C(^0 an d C(^0; the slit f° r C (respectively () variable 
lies along the negative half-axis (if a is real; it is slightly rotated away from the negative 
half-axis if a is complex). However, the orientation is changed (see Figure 3). 

With the above notation for uj and (, the parametrix in U matching N(X) to leading 
order at dll is given by the following expression: 

x 

(69) 



-i(7T/4)(73 V]L ( 1 1 | 7-<7 3 /4 

2 3/2 I -i 1 IS 




(Ti/(A) -, 3 n/2 j ^ 



1 

1 



where y ( = e i7r//2 no;(A). 

Following the steepest descent method, we now formulate a RH-problem for the function 



f $(A)iy(A)-\ AeC\(C/U[/U (-1, 1)) 
i*(A) = < $(A)P(A)- 1 , A e U \ (1 - S, 1] (70) 
U(A)P(A)-\ AeC/\[-l,-l + 5). 

By construction, the function -R(A) has no jumps across (1 — 5, 1) U (—1, —1 + 5). Moreover, 
since apriori R(X) can have no stronger than logarithmic singularities at the points ±1, the 
function R(X) is in fact analytic in the union of the discs U U U. It solves the following 
RH-problem on the contour S = dU U dU U ( — 1 + 5, 1 — 5) (see Figure 4): 

• R(X) is holomorphic for all A ^ S 

• J2(oo) = I 
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Figure 4: Contour E for the i?-RH problem. 



• R+(X) = i?_(A)A(A), A G = £ \ {1 - 5, -1 + 5}, where 

A(x) = N+(x)(l ~ 2f +( x) ^N + (x)-\ xe (-1 + 6,1 -5), (71) 

A(A) =P(\)N(\)-\ XedU\{l-5}, (72) 
A(A) = P(X)N(X)~ 1 , XedU\{-l + 5}, (73) 

Observe that for all — 1 < x < 1 and all < a < 7r, we have 

< /+(x) < e -^=^inf_ 

Moreover, for sufficiently small e there exists a positive constant C5, depending on 5 only, 
such that 

\f+(x)\<e-°; 

for all — 1 + <5 < x < 1 — 5 and all a G V e (7i). Combining the two estimates above, we 
conclude that the jump matrix on [— 1 + 5, 1 — 5] is of order 

I + 0(exp(-C s , E p)), 

where Cs }E is a positive constant which only depends on 5 and e. This estimate is uniform in 

x £ [—1+5,1 — 5], a G [0,7r-e] UX? £ (tt). 

Using (J61|) . we obtain the following asymptotic expansion in inverse powers of ^ I /C ^ for 
the jump matrix on dU: 

P(A)A^(A)- 1 = / + Ai + A 2 + A r , 

A __^(3P 2 -P~ 2 ^W 2 + p- 2 )\ A - JL ( 1 (7 a) 

1 ~ i6vc U 3 /? 2 + r 2 ) -(s/? 2 - r 2 ) y ' 2 ~ 2?c U 1 / ' 

A G dU, 

where /3(A) is defined in (jHEJ)- Since the matrix functions N(X) and A^ _1 (A) are uniformly 
bounded on dU, we conclude from (fKTJ) and (|H2^1 that the error term A r (A) in (J73)) satisfies 
the uniform estimate, 

|A r (A)| < j^j- 2 , (75) 
18 



VA e dU, Va e 



2s 



n 



IT — £ 



UV £ (n), — <So <n. 



Here C$ is a positive constant depending on 5 only. The jump matrix on dU is given by the 
similar representation with the matrices A defined as follows: 



Ai = - 

16 V \ 

A G dU. 



-(3(3- 2 -(3 2 ) i{3f3~ 2 + f3 2 ) 
i(3p~ 2 + /3 2 ) 3(3~ 2 - I3 2 



A.; 



2 7 C 



1 Ai 
-Ai 1 



(76) 



Let us summarize the above calculation. 

Proposition 1 The jump matrix A(A) of the R-RH problem possesses the following proper- 
ties: 



1. For sufficiently small 5 and e, the function A satisfies the estimates: 

|/ _ A( A )| < 9i A G (dU \ {1 — 5}) U {dU \ {-1 + 5}), 
P 



and 



| J - A(x)\ < C s , £ exp(-C s , £ p), x G (-1 + 6, 1 - 5), 
'2s 



(77) 
(78) 



Va G 



7T — £ 



n 



UV e (7c), s <n. 



Here p = n\ sin(a/2)|, and C§, Cs :£ , and C§ )E are positive constants depending on the 
indicated quantities only. The number s is any positive number satisfying the inequality 
Sq > 2/y/5. Moreover, 

P > (79) 



Va G 



2s 



n 



n — e 



2 ' 

UX> E (7r), s <n. 



2. The estimate |77| ) can be extended to the asymptotic series 
fc-i 

A(A) = / + ^A j (A)+A«(A), A G (dU \ {1 — S}) U (dU \ {—1 + 6}), (80) 

where the terms Aj of expansion Il8(ty and the error term Af (A) satisfy the uniform 
estimates: 



|Ai(A)|< 



a 



U) 



I A« (A) I < 



<fe) 



5i) 
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2sr 



n 



71 



UV £ (ti), s <n. 



The positive constants Cg , j = l,...,k are the functions of 5 only. The first two 
terms of the expansion A8U\) . i. e. the functions Ai and A 2 are given by equations \1J$ 
if\EdU\{l- 8}, and by equations $7§ if A G dU \ {-1 + 8}. 

3. Let A u , Arf , and A; denote the limits of A(A) as A approaches the node point 1 — 8 from 
the above, from the below, and from the left along S, respectively. Then these limits 
exist, and the following cyclic equation holds: 



AdA/A„ 1 = /. 

A similar relation (with A; replaced by A r ) holds at the node point —1+8. 



(82) 



4- The matrix function A(A) admits an analytic continuation into a neighborhood of any 
point of the interval (—1 + 8,1 — 8). Moreover, this analytic continuation preserves the 
estimate fl7&| ) with constants Cs )£ , Cs )£ possibly somewhat modified. 



The only statements which need comments are the statements # 3 and #4. These 
statements follow directly from the explicit formulae (j71l - l?B*|) for the jump matrix A(A). 



Corollary 1 The following inequalities hold: 



I/- M 



L 2 (E)nL°°(S) 



< 



k-1 

l J + E A ^- A H i2 



a 



(E)ni«(E) < — f", 



(83) 
(84) 



Va G 



2s r 



n 



7T — 5 



UP £ (vr), s <n, 



where we set Aj = /or A G (—1 + 5, 1 — 8). 



By standard arguments of the L 2 RH theory (see e.g. 0122]), the inequality implies 
the solvability of the i?-RH problem for sufficiently large Sq. Moreover, let Qk, k = 1,2,3 
denote the connected components of the set C\H Then, due to the cyclic relation (J%2~j) . the 
restriction i?|n fc (A) is continuous in Q k for each k (see, e.g., 

To complete the proof of the theorem, we need to show that the solution R(X) of the 
-R-RH problem satisfies the estimates indicated in (|38|). 
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Lemma 1 For sufficiently small 5 and e, and for every k, the function R(X) admits the 
asymptotic representation, 



k-l 



i?(A)=/ + ^it! i (A)+itf)(A), 

3=1 



(85) 



where 



p = n 



a 

sin — 

2 



oo, 



i6) 



uniformly for all A G and a G [^-, tt — el U T> £ [tx). As in the Remark 1 to Theorem 2, 
the latter statement means that there exist positive constants C and Sq such that 



\RjW\ < 77/ 



c 



2s 



l^ fc) (A)| < 



C 



VAGfi (25) , Va6 — .77--; 

The functions Rj{\) are constructed by induction as follows: 

ds 



a + iAi)p fc ' 

UP e (7r), Vn : s < n. 



17) 



RiW = I 



dUUdU S — A 



^ Jauudu s ~ A 



flfc-i(A) = 7^ / E^-i-i.-C^Ai^)-^-, #o s J. (89) 

Z7 ™ JdUUdU j=1 S - A 



Remark 4 We also assume (cf. Remark 1 to Theorem 2) that e is small enough so that 
p > s /2 for all a G [2a, vr - e] U D £ (tt) and n > s . 

The proof of the lemma is essentially a combination of the arguments from [21] and |27j . 

We consider in detail the case of k = 3, which is all that is needed below, but the argument 
extends in an obvious way to any k = 1, 2, . . .. The details are left to the interested reader. 

Write the jump condition for R(X) in the form 



#0+ + Rl+ + R2+ + Rr+ = (R0-+R1- + R 2 - + Rr-)(I + Ai + A 2 + Ar 



(90) 



Here Ai, A 2 are given by (|74|7fiJ) on dU, dU , and we set Ai = A 2 = on (—1 + 5,1 — 5). 
Thus A r = 0(1/ p 3 ) on dU U dU (this error term arises from the Bessel asymptotics: see 
(JHIJ), and A r = 0(e s ' eP ) on (—1 + 5, 1 — 5). We now show that we can define R\ and i? 2 
so that they are of order 1/p and 1/p 2 , respectively. We then show that the remainder R r is 
of order 1/p 3 . Set 

Ro = I. 
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We define Rj by collecting in (|9U|) the terms that we want to be of the same order. First, 

i2i+(A) = i2i_(A)+Ai(A), AGE. (91) 

We are looking for a function i?i(A), which is holomorphic outside E, satisfying Ri(X) = 
0(1/ X), A — > oo, and the above jump condition. The solution to this RH-problem is given 
by the Sokhotsky-Plemelj formula, 

R(X) = C(Ai), (92) 

where 



C(f) = ±-.f /(* 



ds 



s-X 

is the Cauchy operator on E. The condition Ai(A) = 0(1/ p), X G E, p — ► oo (uniform in a), 
implies that there exist c, do, so > such that 

|fli(A)|< (1 + my n ^ s ^ ( 93 ) 

uniformly in a G J^^ 1 , 7r — e] UT> £ (7r) and A satisfying dist(A, E) > d . Actually, this estimate 
is uniform for all A G C \ E up to E. This can be shown either by direct calculation (see 
below) or by suitably deforming the contour E. Indeed, since 

RiW = I AiOO^T, (94) 

JdUUdU S — A 

the estimate (}93|) holds for A up to the interval (—1 + 5', 1 — 5'), for any 5' > 5. Since Ai(A) 
is analytic in the neighborhood of dll U dll (as, actually, are Aj(A) for all j), the contour of 
integration in can be deformed so that the estimate holds up to dll U dll as well. It also 
should be observed that, by the same deformation of the contour of integration in (|94j) . one 
obtains analytic continuations of both the functions Ri + (X) and Ri-(X) in the neighborhood 
of the contour dll U dll. Moreover, this analytic continuation preserves the estimate (JOB}. 

Now define ^(A) by the jump condition 

i? 2 + (A) =i? 2 ^(A) + J Ri_(A)Ai(A) + A 2 (A), A G E, (95) 

together with the requirement of analyticity for A G C \ E, and the condition i^(A) = o(l) 
for A — > oo. The solution to this RH-problem is 

R 2 (X) = C(R 1 ^A 1 +A 2 ) = 

I (R 1 -(s)A 1 (s)+A 2 (s))^-, XeC\(dUUdU). (%) 

^ l JdUUdU S — A 

Using ()93p. the estimates Aj = 0(1/ fP), and the analyticity of and Aj in the neighbor- 
hood of dll U dll, we obtain in the same way as for R±. for some c > 

WA)|< (i + | A |)p2 » AGC\E, n>s (97) 
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with the same uniformity and analyticity properties in a and A. Below in the proof, the 
same symbol c will stand for various constants independent of a, A, and n. 

Now from ()90l91l95j1 we obtain 

Rr + (X) = M(A) + i? r _(A)A(A), AGS, (98) 

where 

M = R 2 ^A 1 + _ + i? 2 _)A 2 + (I + Rx _ + R 2 -)A r . 

Remark In the terminology of [21], equation ()98j) is an inhomogeneous RH-problem of 
type 2. 

Since R r = R — I — R± — R 2 , the matrix function R r (X) is holomorphic outside S and 
satisfies the condition R r (X) = o(l) as A — > oo. Therefore, 

R r (\) = C{M) + C(R r -(A-I)), AgC\S. (99) 

(It is worth mentioning that, by virtue of property # 3 of the jump matrix A(A) formulated in 
proposition ^ equation (|9~9~|) is consistent with the absence of the singularities of the function 
R r (X) at the node points 1 — 5 and —1 + 5.) Equation (|9*9*|) . in turn, implies that 

R r _(\) = C-(M) + C-(R r -(A-I)), AGS, (100) 

where C_(/) = limy^,\C(f), as A' approaches a point A G S from the "— " side of S. Now 
defining the operator 

C A (/) = CL(/(A-J)), 

we represent fllOOj) in the form 

{I - C A )(R r _) = C_(M). (101) 



Because of the L°° part of the estimate ()84|) . and the fact that C_ is a bounded operator 
from £ 2 (S) to L 2 (S), it follows that the operator norm ||Ca|| = 0(1/ p), and hence / — C\ 
is invertible by Neumann series for sq (and therefore p) sufficiently large. Thus (jlOlj) gives 



i? r _ = (/-C A )- 1 (C_(M)). (102) 

UL,L( ' "i«.< in -i ■" /||l-i>:.i ■-- "' 

Together with ()102|) . this yields the uniform estimate, 



Moreover, using the L 2 part of the estimate we conclude that ||C_(M)|| i 2( S ) = 0(p 



Rt-\\l^)<^, (103) 



Va G 



2s 
n 



UV £ (tt), n>s . 



Combining the estimate ()103j) with equation (j9Uj) . we can complete the proof of the lemma 
as follows. 
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First, assuming that dist(A, E) > do, we immediately arrive at the estimate 

|g(M)(A)l< ^ n>s , (104) 

for the first term in the r.h.s. of ()99|) . and the estimate 

|C( J R r _(A-/))(A)| < ^-f—lli^.ll^llA-Jll^s) < 

c 



n > s , (105) 

)th the estimates are uni 

they yield the estimate 



(1 + |A|)P 4; 

for the second term. Both the estimates are uniform in a G [— , n — e\ U V £ (n). Together 



|flr(A)|< , n>s , (106) 

uniformly in a G [— , vr — el U V e (ir) and A satisfying dist(A, S) > do. 

Second, we observe that the matrix A r (A) coincides with the matrix A(A) — / on the 
interval (—1+5, 1—5). Hence, by property # 4 of the matrix function A(A) (see proposition^), 
the matrix function A r (A) admits an analytic continuation in the neighborhood of any point 
of the interval (—1 + 5,1 — 5), and this continuation preserves the estimate, A r = 0(e~ Cs - eP ). 
This means that, by bending the segment (— 1 + 5, 1 — 5) of the contour E we can extend A 
in the estimate ()104j) up to the interval (—1 + 25,1 — 25). Using property # 4 of the jump 
matrix A(A) one more time, we can rewrite the second term in the r.h.s. of equation (|99j) as 



C(J2 P _(A-J))(A) = -^/ R r 4s)(A(s)-I) <h 



A 



If ds 
R r (s)(A(s) - 1)- 



(107) 



2vri 7 7 (d) s - A 



if A lies above the interval (— 1 + 25, 1 — 25), and as 

c(iMA - *))(A) = ^- [ Rr-(*)(M*) - J )- 

zm JdUUdU 

-L/ (R r (s) - M(s))(I - A-\s))- 
2m J ( U ) 



ds 



A 



ds 



:io8) 



if A lies below the interval (—1 + 25,1 — 25). Here, the contours 7^ and 7^"- ) are the 
slight deformations of the segment (—1 +5,1—5) down and up, respectively. Using, in 
representations ()107|) and ()108j) . the estimate ()10ti|) for R r (X), we extend the variable A in 
the estimate (|l()5j) up to the interval (—1 + 25,1 — 25). 

The above extensions of the estimates ()104j) and (J105)) mean, in particular, that they 
both, and hence the estimate f)lU6j) . are valid for all A G ^ 25 \ The proof of the lemma is 
completed. □ 
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We now derive explicit formulae for the terms -Ri(A) and i?2(A) of the expansion ()85|). 
By Lemma 1, 



m) = J_[ Me>£ r 2W = ±[ fi-MAiM + A^ (109) 

As noted in [27|, we can also obtain the expressions for Rj{\) in the following way. It is 
not difficult to check that Ai(A) and Az(\) are analytic in (U \ {1}) U (U \ { — 1}) with the 
simple poles at ±1. We have 

AO) BO) 
Ai(A) = t — - + 0(1), asA->l, Ai(A) = — — + 0(1), as A -+ -1, (110) 

A — 1 A + 1 

where the constant matrices AO) and BO) are obtained by expanding ui(\) and /5(A) in 
(I74l76j) at A = ±1. It is easy to verify directly that the Riemann-Hilbert problem for R\(\) 
has the solution: 



Rl (X) = < ^ + ^ > iorXeC\(UUU) 
l[) ^S + fS-A^A), forXeUUU. [ j 

Using the series ()54|) and the expansion of /3(A) at ±1, it is not difficult to obtain the 
singular and constant term in the Laurent expansion of Ax (A). By the first formula in (Jlllj) . 
we obtain (using the singular term) the expression (I39j) . 

Similarly we may calculate the singular term in the expansion of A2(A) at ±1, and use 
the second formula in (jlllj) to evaluate i?i(±l) (note that the formula (j39j) is valid only 
outside UUU). It is then easy to compute the integral for R 2 in (j!09j) and obtain (j4"Uj) . This 
completes the proof of the theorem. □. 

Now we give some remarks and corollaries of Theorem |21 
Remark 5 Estimate (j37j) and formula (j36j) imply that 

0(n, a) ~ cos — , A(n, a) ~ sin 2 — , 

and using either (j3*3*j) or (j3*4"j) we recover the master term of Widom's asymptotics p a (cf. 
also |7j), 

n i a 

mD n \a)~n In cos— , n — >• oo. 



Corollary 2 The function A(n, a) admits the asymptotic expansion 

a / \ . 9 ol cos 2 (a/2) / 1 \ . 9 . 

A(n,a) = sm 2 J_L_i + o ( — sin 2 a, p -> oo, (112) 

2 4n 2 VP/ 
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which is uniform for a E [^-, 7r] . 
Remark 6 The statement, 



A r (n, a) = O I — ) sin a 

P 3 



p — > oo 



uniformly for a G [— , 7r] , means that there exist positive constants C and Sq, such that 

C 



\A r (n, of) < — r sin a, 



113) 



Va G 



2s r 



7? 



7T 



and So < n. 



Proof of Corollary 2. To calculate A we need the asymptotics of $(A) outside the neigh- 
borhoods U and U . By ()38|) these are given by the expression: 



§{\) = {I + R l + R 2 + Rf ) )N, \ett {25) 



where Rf 1 is estimated by (jS7j) for k = 3. In particular, the estimate (J87)l becomes 



a 



a 



0(p Ism — if A = ±zcot— . 



Similarly, 



Since 



a 



Rj(±i cot a/2) = 0(p~ 3 ) sin -, j = 1, 2, . . . 



iV(±i cot a/2) 



cos(a/4) ±sin(a/4) 
=Fsin(a/4) cos(a/4) 



(114) 

(115) 

(116) 
(117) 



we have N(-i cot a/2)' 1 = N(icota/2) = 0(1) and [N(-i cot a/2)~ l N(i cot a/2)] 12 
sin a/2. Definition (jHSJ) and equations (|114J115Jllfi|) then imply 



A(n, a) = [iV(-i cot a/2)" 1 (l + (O(p^) + 0(p- 2 ) + 0(p~ 3 )) sin |) 



/ + (CKp- 1 ) + 0(p- 2 ) + 0(p~ 3 )) sin | J iV(i cot a/2) 
= A (a) + -Ai(a) + -^A 2 (a) + 



12 



where 



a 



— sin (a/2), 



A (a) = sin — , 



118) 



(119) 



(120) 



and f(a, n) is uniformly bounded for a G [2s /n, n — e] U P e (7r), and s < n. Note that to 
write ()118|) we used the fact that det R(X) = 1. 
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In order to determine the terms Ai(a) and A 2 (a) in this equation we need Ri^(±i cot a/2). 
These values we obtain from (|39|4O0 : 

fl l( ±ico ta /2) = ±± ( - C °f{*> , (121) 

8n \±sm(a/2) cos(a/2) J 

R 2 (±i cot a/2) = ±-L ( ± - 8 COS / a / o 2 ) "I . (122) 

v 7 ; 2 7 np V 8cos(a/2) ±sin(a/2) J v ; 

As det$(A) = 1, the inverse $ _1 is easy to compute, and after a simple computation we 
arrive at the equations, 

Ai(ct) = 0, A 2 (a) = -- cos 2 (a/2). (123) 

Now it only remains to show that the function f(a,n) in ()119|) satisfies the estimate 

\f(a,n)\ < Ccos 2 (a/2), (124) 

for a G [2s /n, n], and s < n. In fact, since the uniform boundedness of f(a, n) on the sets 
indicated has already been established, it is enough to show that estimate ()124|) holds for all 
a G T> e /4(ir). 

Observe that, for fixed n, the quantity A(n, a) is an analytic function at a = it with 
A(n, 7r) = 1 and (d/da)A(n, it) = (see ([170 and (jSJ)) so that we can write down the Taylor 
series for A in j3 = it — a: 

A(a) = 1 + (a - n) 2 a 2 + ■■■■ (125) 

As follows from equation ()119|) . f(a,n) is a holomorphic function of a in T> £ (tt). Using a 
representation of f(a, n) by a Cauchy integral, we obtain: 

f(a,n) = — — da = — — da+ 

2m J&D e/2 (n) " - " 2 ™ Jdv, /2 (n) " - 

Q . / 7^ T^a+ . / r^c rda, 7r-a <e/4. 

(126) 

At the same time, from ([1190 . ([1250 . and ([1230 it follows that the Taylor series of f(a, n) at 
a = Ti has the form, 

f(a, n) = (a — 7i) 2 d 2 + . . . 

Therefore, the first two integrals in the r.h.s. of ([1260 must be zero, and the third one, by 
virtue of the uniform boundedness of f(a,n) for all a G T> £ (jr) and all n > sq, yields the 
estimate ([1240 for all a G V e /±{it) and all n > sq. The proof of the corollary is completed. 
□ 



Remark 7 Here is an alternative derivation of the leading terms in formula ([1120 . 
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We start with equations (jll9ll20j) . The issue is the exact evaluation of the quantities 
Ai(a) and A 2 (a). This can be done with the help of the relevant (integrable) differential 
system associated in the standard way with the original m - RH problem. Indeed, it is shown 
in that the Toeplitz determinant D n (a), considered as the function of the variable 



t = e~ 2ia 



is the t - function for the Painleve VI equation characterized by the parameters 

#00 = —9o — n , 9i=9 t = 0, 



where we use the ^-notations of Jimbo, see [28J. According to [28], this means that the 
quantity 

V (t)=t(t-l)j t \nD n (127) 
satisfies the following nonlinear differential equation (the r- form of Painleve VI): 



+ 



dri n 2 \ ( dri \ ( dri\ n 2 dri 

t— -71 - — + 



2 . . 4 



dt 4 J \ dt J \dt J 2 dt 
The functions A(n,a) and rj(t) = r](n,t) are related by the equation 



drj 
~dt 



(128) 



^ 1 — t drj ^ 1 (129) 
n 2 dt n 2 ' 

and we may anticipate an expansion for 77 similar to ()112|) . Indeed we expect 

r l (t) = ri(n,t)=n 2 r ]0 (t)+nri 1 (t)+r ]2 (t) + o(^ S j , p -> 00, (130) 

where 

Vo (t) = ^(l-Vt) 2 . (131) 

A substitution of the asymptotics ()13()j) into the equation ()128|) gives us, after a straightfor- 
ward calculation, the following formulae for the coefficient functions t]i(t) and 7]2{t): 

rh(t) = 0, m(t) = -^(l + Vi) 2 . (132) 
These equations together with p29j) lead immediately to the leading terms in the formula 

dm. 

It also should be noticed that the differentiability of the asymptotics ()119j) follow from 
its uniformity in the disk V e (7i). 
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4 Asymptotic evaluation of D n (a). Proof of estimate 

OH). 



The asymptotic evaluation of the Toeplitz determinant D n (a) is based on the integration of 
the differential identity (|34|) from a to ao (which is close to tt from below). We have: 



(a - a)(\nD n )'(a ) - \nD n (a ) + \nD n (a) = - 




(133) 



Fix n and set ao = vr — (3. Substituting for lnD n (ir — (3) the expansion (jl7j) . and for 
A(0) the asymptotics pi2jl . and after taking the limit (3 — > 0, we immediately obtain (fT2*|) 
with the remainder 0(l/{nsin(a/2)}) uniformly for — < a < n — e, n > so, e > 0. 
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